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Abstract 

In the framework of the variational principle there are introduced canonical variables 
describing magnetohydrodynamic (MHD) flows of general type without any restrictions 
for invariants of the motion. It is shown that the velocity representation of the Clebsch 
type introduced by means of the variational principle with constraints is equivalent to 
the representation following from the generalization of the Weber transformation for the 
case of arbitrary MHD flows. The integrals of motion and local invariants for MHD are 
under examination. It is proved that there exists generalization of the Ertel invariant. It 
is expressed in terms of generalized vorticity field (discussed earlier by Vladimirov and 
Moffatt (V. A. Vladimirov, H. K. Moffatt, J. Fl. Mech., 283, pp. 125-139, 1995) for the 
incompressible case) . The generalized vorticity presents the frozen-in field for the barotropic 
and isentropic flows and therefore for these flows there exists generalized helicity invariant. 
This result generalizes one obtained by Vladimirov and Moffatt in the cited work for the 
incompressible fluid. It is shown that to each invariant of the conventional hydrodynamics 
corresponds MHD invariant and therefore our approach allows correct limit transition to 
the conventional hydrodynamic case. The additional advantage of the approach proposed 
enables one to deal with discontinuous flows, including all types of possible breaks. 

1 Introduction. 



It is well-known that description of the solid media flows in terms of the canonical (hamil- 
tonian) variables is very useful and effective, see for instance fl|, yj. In terms of the hamiltonian 
variables it is possible to deal with all nonlinear processes in unified terms not depending on the 
specific problem related to the media under investigation. For instance, all variants of the per- 
turbation theory are expressed in terms of different order nonlinear vertices, which along with 
the linear dispersion relation contain the specific information relating to the concrete system 
under investigation, cf. Refs. [0, In the problems of the nonlinear stability investigations the 
conventional Hamiltonian approach based upon the corresponding variational principle allows 
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one to use the Hamiltonian along with other integrals of motion (momentum, number of quasi- 
particles, topological invariants) in order to construct the relevant Lyapunov functional, cf. Refs. 
H 0) H) Hi- Therefore, it makes important the problem of introducing the canonical variables 
and corresponding variational principle for the general type MHD flows (i. e., non-barotropic 
and including all types of the breaks possible for MHD) and obtaining the complete set of the 
local invariants, see definition and discussions in original papers (T^, ITT], [TJ], [i3|] and in the recent 
review [|TJ . As for the first item, the example of the variational principle describing all possible 
breaks is presented in the recent work jl4| . 

Here in the framework of some modification of the variational principle of the cited work we 
examine the problem of the MHD invariants. Note that the set of invariants for MHD discussed 
in the literature till now is incomplete. It becomes evident if one takes into account that for 
the vanishing magnetic field this set has to go over to the set of the conventional hydrodynamic 
invariants. But this limit transition does not reproduce Ertel, vorticity and helicity invariants 
existing for the hydrodynamic flows. For the particular case of incompressible MHD flows 
generalized vorticity and helicity invariants were obtained in the paper ||. Below we show that 
the generalized vorticity and helicity invariants exist also for compressible barotropic flows, and 
derive MHD generalization for the Ertel invariant. 

The plan of the paper is as follows. In section |2] we briefly discuss appropriate variational 
principle, introducing the Clebsch type velocity representation by means of constraints and 
defining the canonical variables. In the following section |3| we develop generalization of the Weber 
transformation and show that it leads to the velocity representation, which is equivalent to the 
one following from the variational principle under discussion. In section |] we examine MHD 
integrals of motion, introducing 'missing' MHD invariants, and discuss their transformation 
properties relating to change of gauge. In section [5] we make some conclusions and formulate 
problems to be solved later. 



2 Variational principle and canonical variables. 



Let us briefly describe the variational principle and subsidiary variables describing dissipation- 
free MHD. Starting with the standard Lagrangian density 

v 2 H 2 

where p, s and e(p, s) present the fluid density, entropy and internal energy, respectively, H 
denotes magnetic field, we have to include to the action A the constraint terms. Then the action 
can be presented in the form 

A = J dW , L' = J drC , £' = £ + C c , (2) 

where £ c is the part of the Lagrangian density respective for constraints, 

/ dA \ HcurlA 
£ c = pDcp + \Dn + aDs + M • — - [v, curl A] + VA + — . (3) 

\ Ot ' J 4:71 
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Here D = d t + (vV) is substantial (material) derivative and A is the vector potential.^ Including 
the last two terms into C c allows us to introduce relation H = curl A strictly into the variational 
principle (it follows after variation with respect to H). 

Supposing first that all variables introduced (including velocity) are independent we obtain 
the set of variational equations of the form 

5<p =>■ dtp + div{p\) = 0, (4) 



§p ==> Dcp = w- v 2 /2 , (5) 

5X^ Dfi = 0, (6) 

5p m d t \ m + div(\ m v) = 0, (7) 

5(7 =7* Ds = 0, (8) 

5s =>• d t o + div(av) = —pT, (9) 

5M d t A = [v, curl A] - VA, (10) 

5A^ d t M = ^-^ + curl[v,M]. (11) 
An 

5H^ H = curlA, (12) 



SA ==► dzwM = 0, (13) 

where u; and T are the enthalpy density and temperature. 

Note that in this section we suppose the velocity field to be independent on other variables. 
Therefore, variation with respect to v gives us the velocity representation: 

5v pv = -pVip - X m Vp m - crVs - [H, M] . (14) 

It is convenient to rewrite it in a shortened form that emphasizes it's structure. Bearing in mind 
that the velocity potential ip, vector Lagrange markers /la, entropy s and the vector potential A 
can be treated as generalized coordinates, one can see that p, A, o and subsidiary field M are 
conjugated momenta, respectively. Let 

Q=(Q,A), Q = (v?,m,s), V = 5A/5d t Q, V = (P,M)- (15) 

1 This form of the action slightly differs from the one proposed in Ref . jl4| . The main difference consists 
in introducing the vector potential for the magnetic field. Therefore, here the canonical pair is A, M instead of 
H, S, where S = citr/M. We do not deal here with the discontinuous flows and thus we omit the surface term in 
the action. But adding corresponding surface term we can easily take the breaks into account. 
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Then the velocity representation takes the transparent form 



v = v (?,VS), Vo = v ft + v M , v h = --VQ, v M = --[H,M]. (16) 

P P 

Here subindexes h and M correspond to the "hydrodynamic" and "magnetic" parts of the ve- 
locity field. The hydrodynamic part, v^, corresponds to the generalized Clebsch representation, 
compare papers [113], P] , and magnetic part, Vm, coincides with the conventional term if we 
replace the divergence-free field M by curlS. The latter was first introduced by Zakharov and 



Kuznetsov, cf. Ref. 17 



From the velocity representation Eq. fllED and the equations of motion (|4j) - ([11]) it strictly 
follows that the velocity field v = v satisfies Euler equation with the magnetic force taken into 
account. Namely, providing differentiation we obtain 

p Dv = -Vp + — [curlH, H], (17) 

where p is the fluid pressure. 
2.1 Canonical variables. 

The variational principle can be easily reformulated in the Hamiltonian form. Excluding the 
magnetic and velocity fields by means of Eqs. (|T2|), flT6|) we arrive to the following Hamiltonian 
density 

H = H(V, VQ) =Vd t Q-C' = / & + pe(p, s) + ^£ - (M, VA). (18) 



Equations of motion (|4]) - ([II]) can be expressed now in the canonical form 

d t Q = 5H v /6V, d t V = -5Hv/5Q, Q = (y?, //, s; A) , V = (p, A, a; M) ; (19) 



Eq. ( p!2|) serves as a definition of magnetic field, and divergence-free condition for the subsidiary 
field M, Eq. (|l^), follows from variation of the action 

A = J dt J dr (Vd t Q - H) (20) 

with respect to A. Note that it is possible to put A = 0. Under this assumption the divergence- 
free condition for the field M vanishes, but from Eq. ([11]) follows that divM. is conserved quantity, 
and supposing that &M = holds for some initial moment we arrive to the conclusion that 
this is valid for arbitrary moment. Nevertheless, it proves convenient to deal with A ^ that 
makes it possible to use different gauge conditions for the vector potential. 

The variational principle presented gives us the set of dynamic equations from which follow 
conventional MHD equations, (|4j), (|S]), flT?P and equation for the frozen-in magnetic field, which 
follows from Eq. fllCf) after taking curl operation, 

d t U = curl[v,B]. (21) 
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On the contrary, if at some initial moment, t = t , we have the conventional MHD fields p, 
s, v and H, then we can find the initial subsidiary fields <p, fit, A, a, A, M and A, satisfying 
Eqs. (|T2"D - (|14"1) . This can be done up to the gauge transformations (do not changing the velocity 
and magnetic field) due to the fact that the subsidiary fields play a role of generalized potentials. 
Then, if the uniqueness conditions hold both for the conventional MHD equations and for the 
set of variational equations, we arrive to the conclusion that corresponding solutions coincide for 
all moments. In this sense we can say that these sets of equations are equivalent, cf. Ref. ||. 

The complete representation of the velocity field in the form of the generalized Clebsch 
representation, Eq. (|T6| ) allows, first, to deal with the MHD flows of general type, including all 
types of breaks, cf . Ref. |14| ; second, for the zero magnetic field it gives correct limit transition 



to the conventional hydrodynamics, cf. Refs. |nj, |15 |; third, it allows to obtain the additional 



to the known ones integrals and invariants of motion for the MHD flows: for instance, generalized 
Ertel invariant, generalized vorticity and generalized helicity, see below. The two last integrals 
were deduced for the particular case of incompressible flows in the paper ||. 

Moreover, it is possible to show that representation ([16]) is equivalent to the one following 
from the Weber transformation, cf . Refs. [18|, [19|] and the recent review U . The generalization 
of the Weber transformation for the ideal MHD incompressible flows was obtained by Vladimirov 
and Moffatt, cf. Ref. B. 



3 Generalized Weber transformation. 



Suppose here that the fluid particles are labelled by Lagrange markers a = (ai,a 2 ,a 3 ). The 
label of the particle passing through point r = (x^a^,^) at time t is then 

a = a(r,t), (22) 

and 

Da=^ + (vV)a = 0. (23) 
ot 

The particle paths and velocities are given by the inverse function 



<9r 

r(a,t), v = Dr(a,t) = — 



(24) 



a=const 



Let the initial position of the particle labelled a is X, i.e., 

r(a,0)=X(a). (25) 

A natural choice of label would be X(a) = a; however it is convenient to retain the extra freedom 
represented by the "rearrangement function" X(a). 

We now seek to transform the equation of motion flT7| ) to integrable form, by generalization 
of the argument of Weber |18[ (see, for example, Refs. [^(J, []]], and ||. It is convenient to 
represent here the equation of motion in the following form 

Dv = -Vw + rVs + [J,h], (26) 
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where h = H/p and the vector J is defined according to 

curl~H. 



(27) 



(Dv K ) dXk = 9W 9Xk I T 9s 9Xh I [J h] K dXk (28) 
dcii dxk da>i dxk ddi ' ddi 
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being proportional to the current density. Multiplying Eq. (1261) by dxk/ddi we have 



The l.h.s. can be represented in the form 



dxk ^ ( dxk \ d , 2 



^•^r D K ,k ^)-^yi % (29) 

where we have taken into account that operator D = d / dt\ a=const and therefore Dxk = ffc and 
D commute with derivative 8/ ddi. Eq. takes now the form 

dxk \ d , 9 . . m <9s r ., , , dxk 



It is convenient to transform the last term by means of the dynamical equation for the subsidiary 
field m = M/p (compare Eq. (|TT|)) 

Dm = (m, V)v + J /p. (31) 



Then we can transform the last term in the r.h.s. of Eq. fl3"0|) to the form of substantial derivative, 
see Appendix 

lJM k ^ = D([ m ,H} k ^). (32) 



Odi \ Odi r 

Analogously, the first two terms in the r.h.s. of Eq. ([!(]) can be presented as substantial 
derivatives by means of introducing subsidiary functions <p and a, which satisfy equations (com- 
pare Eqs. ©, ©) 

D (jj = —T, (33) 
Dip = w - v 2 /2. (34) 

Then 

^ ddi da^ \pj ^ \ddi p J J dd^ V ^ W ^ ^ \ddi) ' 

where we have taken into account that Ds = along with D(ds/dd,j) = 0. Therefore, we can 
present the Euler equation fl30|) in the integrable form 



D I v k p) = - D (¥)-d(£--)+d( [m,H]^ ) ■ (36) 



ddi J \ddi J \ddi p J \ 3d 
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Integration leads to the relation 

dx k <9<p ds a dx k 

Vk^~ = --5 ^ H > m fc a h b ^ ( 37 ) 

da, L oai oai p oai 

Here b = b(a) does not depend on time explicitly, Dh = 0, presenting the vector constant 
of integration. Multiplying this relation by dai/dxj allows reverting from Lagrangian variables 
(a, t), to the Eulerian ones, (r, t), 

v = -Vf + b k Va k - -Vs - [h, M]. (38) 

P 

This representation obviously coincides with the Clebsch representation obtained above from 
the variational principle with constraints if one identifies b with —A/ p and a with \i. Moreover, 
this proves equivalence of description of the general type magnetohydrodynamic flows in terms of 
canonical variables introduced and the conventional description in Lagrange or Euler variables. 
The equations of motion for the generalized coordinates and momenta follow now from definitions 
of the subsidiary variables a, m = M/p, a, <p and b. 

Emphasize here that the vector field M = pm introduced by Eq. ([31]) satisfies integral relation 

d t [ (M,dE) = / (J,dE), (39) 
is Jt, 

where S presents some oriented area moving with the fluid. The proof of this statement see 
in Appendix. Expressing M = curlS and making use of the Stocks theorem we arrive to a 
conclusion that time derivative of the vector S circulation over the closed frozen-in contour <9£ 
is proportional to the current (remind, J = (47r) _1 cur/H and differs from the current density by 
constant multiplier) intersecting the surface defined by this contour, 



dt [ (S,di) = [ (J.tffl) = (An)- 1 [ (H,dl) 



(40) 



that highlights the physical meaning of the subsidiary field S usually introduced for the canonical 
description of MHD flows. This fact was first indicated in Ref. for the incompressible flows. 
Now we see that it holds true for the general case. 

The vector constant of integration, b, may be expressed in terms of the initial conditions, 



<p Q = <p(a, 0), c 



h T7 ( \ dXk j_ j_ ds 
h = Vfc(a)— h ^ h Co— , 

oai oai oai 
, F*(a) = V fc (a) + [ho,Mo]*, K fe (a) = v*(a, 0) , ^ 

t=o 



h = h (a) = h(x(a, 0), 0) = h(X(a), 0), M = M (a) = M(x(a, 0), 0) = M(X(a), 0). 
Under special conditions, namely, for 

X(a) = a, r(a, 0) = a, a(r, 0) = r, (42) 

from Eq. (f4~l|) follows 

b^V^ + ^ + Co^. (43) 
oai oai 
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Adopting zero initial conditions, 



M = 0, Vo = 0, a = 0, (44) 

we obtain 

b = V(a) = v(a, 0) = v (a) = v(a, 0) , (45) 

where tilde indicates that we are dealing with the velocity field in the Lagrange description, i.e., 
v(a, t) denotes velocity of the fluid particle with label a at time t. Evidently, v(a, t) = v(r, t), 
where a and r are linked by relations fl22|) and (^4|) for the specific choice given by Eqs. ([E|), 
fl44j) . Then the velocity representation takes the particular form 

v = v h -[h,M], w h = -X7cp + v ok Va k -?-X7s, v (r) = v(r, 0), a(r, 0) = r. (46) 

It differs from the one presented in Ref. [|j] by involving the entropy term. Emphasize here 
that existence of this term allows to describe general type MHD and hydrodynamic flows with 



arbitrary possible discontinuities, including shocks, slides and rotational breaks, cf. Ref. [16 
One can omit this term for continuous barotropic and isentropic flows. 



4 Integrals of motion. 



The conservation laws, as it is well-known, follow from the specific symmetries of the action. 



Existence of the relabelling transformations group (first discussed by Salmon in Ref. pl|) of 
the Lagrange markers, fi 1 leads to the additional to the energy, the fluid momentum and mass 
integrals of motion. These additional integrals are expressed in terms of the Lagrange description 
of the motion, i.e., in terms of the Lagrange markers, etc. Therefore, as a rule, they are gauge 
dependent. The frozen- in character of the magnetic field leads to the specific topological integrals 
of motion, namely, magnetic helicity and, cross-helicity, first discussed by Moffatt in Ref. |[22| , 
see also review |3J. Corresponding densities are respectively 

/i M = (A,H), (47) 

and 

fcc = (v,H). (48) 

As it strictly follows from the dynamic equations, the local conservation law for the magnetic 
helicity holds true for general type MHD flows 

d t h M + divciM = 0, qjvr = vh M - H ■ ((A, v) + A) . (49) 

On the contrary, the cross-helicity in general case is governed by equation 

d t hc/dt = -div [vh c + (w- v 2 /2)B] + Tdiv(sH) 

and is not conserved. But for barotropic and isentropic flows the pressure p = p(p) and he is 
conserved: 

dthc + divqc = 0, q c = vh c + (X- v 2 /2)H_, (50) 
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where x — f dp/ p. 

For the general case there is known one more conserved quantity first discovered by Gordin 
and Petviashvili, cf. Ref. f23|. Corresponding density is 



/i P = (H,Vs), (51) 

and 

d t hp + divqp = 0, qp = vhp . (52) 

With this local conserved quantity there is linked integral conservation law. Namely, inte- 
grating hp over arbitrary substantial volume V we obtain conserved quantity Ip, 

J P = drhp , d t Jp = 0. (53) 



v 



Note here that the latter quantity gives us example of the so called local Lagrange invariants, 
cf. Refs. [II], [II, [12], |13[ and @, |I|. By definition they obey the following equations 



d t a + (v, V)a = , d t I + (v, V)I = , (54) 
<9 t J + (v, V) J - (J, V)v = , (55) 

d t L + (v, V)L + (L, V)v + [L, curly] = 0, or, equivalently, d t L + V(v, L) - [v, curlL] = 0. (56) 

Here a and I present the scalar and vector Lagrange invariants, J is frozen-in field, and L 
presents S"-type invariant in terminology of Ref. [fL2l , related to frozen-in surface. To these 
invariants it is necessary to add the density p. Evidently, the quantity hp/p presents a-type 
invariant. The Lagrange markers /J, and quantities X/p give us examples of the vector Lagrange 
invariants, magnetic field H divided by p, h = H/p is invariant of the J- type, gradient of any 
scalar Lagrange invariant is S"-type invariant, 

L' = V«. (57) 

There exist also another relations between different type invariants, see Refs. 0, [I[, allowing 
to produce new invariants. For instance, scalar product of the J and L invariants presents some 
scalar Lagrange invariant, symbolically 

a' = (J,L). (58) 

The presented above invariant hp/p can be obtained by means of this relation if we put J = h 
and L = Vs. Another examples present relations generating J- (L)- type invariants by means 
of two L- (J-) type invariants, 

J' = [L,L']/p, (59) 
L' = p[J,J']. (60) 
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Note here that integrating of the density Hm over arbitrary substantial volume does not lead 
to the conserved integral. It is easy to check that 

1m = / drh M (61) 
Jv 

satisfies relation 

d t l M = [ dS((A,v) + A)H n , H n = (H,n), (62) 
Jav 

where integration in the r.h.s. is performed over the boundary dV of the volume V, n is 
outward normal and oE presents infinitesimal area of the surface dV. It is obvious that Xm 
will be integral of motion if H n equals zero. This fact leads to the conclusion that Xu presents 
integral of motion if we choose the substantial volume in such a way that initial volume V|t=t 
is such that H n \ t= t = because this condition is invariant of the motion: if equality H n = 
holds for the initial moment then it holds true in the future. 

Another way to make Xm invariant consists in fixing the gauge of the vector potential A in 
such a way that (A, v) + A = 0. Then the dynamic equation for A, (|1(J), takes the form 

d t A + V(v, A) - [v, curl A] = 0, 

i.e., A presents invariant of the L- type. Under this gauge condition the quantity Hm/p presents 
the scalar Lagrange invariant, D(Iim/p) = 0. 

As for the local conservation law for the cross-helicity, Eq. (|50"D , it obviously leads to the 
integral conserved quantity Xq for the barotropic flows but with following restriction: integration 
have to be performed over the specific substantial volume such one that condition H n \ g y = 
(this condition is invariant of the motion) holds, 

d t X c = , X c = drh c , H n \ d y = 0. 
Jv 

Existence of the recursive procedure allowing one to construct new invariants on the basis 
of the starting set of invariants, see Refs. 0, [![, accentuates the role of the local invariants 
among other conserved quantities. Although in terms of the Lagrangian variables (such as the 
markers fx) there exist a wide set of invariants, see, for instance, Ref. |IJ, the most interesting 
invariants are such that can be expressed in Eulerian (physical) variables and are gauge invariant. 
Emphasize here that in the conventional hydrodynamics there exists Ertel invariant ag, 

a E = h E /p, h E = (<jj,Vs), (63) 

where uj = curlw is vorticity, 

d t h E + divq E = 0, q E = h E v, Dct E = 0. (64) 
Corresponding integral of motion reads 

d t X E = 0, X E = [drh E . (65) 
Jv 

Note here that X E = holds true for arbitrary substantial volume V. 
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The Ertel invariant density has the structure of the Eq. ( j58|) with L = Vs, J = o?/p, where 
uj is vorticity, uj = curlv (remind that a; is a frozen-in field for the barotropic hydrodynamic 
flows). In the hydrodynamic case there exists also the helicity invariant 

h H = (w, v), (66) 

which has topological meaning, defining knottness of the flow. It satisfies equation 

d t h H + divq H = 0, q H = h H v + (x ~ v 2 /2)uj, (67) 

and evidently leads to the corresponding integral conservation law 

d t X H = , for u n \ g y = , X H = drh E . (68) 

Jv 

For the MHD case the vector uj/p does not present frozen-in field due to the fact that 
magnetic force is not potential. It seems rather evident that for the MHD case there have to 
exist integrals of motion generalizing the conventional helicity and Ertel invariant along with 
vorticity one, which have to pass into conventional ones for vanishing magnetic field. The 
generalization for the vorticity and helicity invariants was obtained by authors of the paper || 
for the particular case of the incompressible flows. In the following section it is shown that there 
exists MHD generalization for the Ertel invariant, and results of the paper || relating to the 
vorticity and helicity can be extended for incompressible barotropic MHD flows. 

4.1 Generalized vorticity. 



Let us prove that the quantity u: h /p, where 



uj h = curl\ h = - 



V(^],VQ 



V ( 

P 



V/i r , 



p 



(69) 



presents frozen-in field ('hydrodynamic' part of the vorticity) for the barotropic MHD flows. It 
would be trivial consequence of the fact that [L, L']/p, where L, L' are Lamb type invariants, is 
local invariant of the frozen-in type if all quantities Q and Pj p satisfy homogeneous transport 
equations being a- or I type invariants (remember, that Va and V/ m are L type invariants). 
But (p and a/p satisfy inhomogeneous equations of motion. Therefore, let us start with equation 
of motion for the 'hydrodynamic' part of the velocity. Differentiating representation (|T6|) and 
making use of relations 

D{VX) = W(DX) - (Vv ro ) • d m X 



we have 



D [ - J • VQ - - • V{DQ) + -(VO • d rn Q = TVs - V(u> - v 2 /2) - v hm (Vv m ) 
P ) P P 



Dv h ~- 

or, after simple rearrangements 



Dv h = -Vp/p + (v m - v hm ) ■ VX 



(70) 
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Taking the curl of this equation leads to 

d t uj h = -curl((v m d m )v h ) + [Vp, Vp}/ p 2 - curl(v hm Vv m ) = 

= [Vp, Vp]/p 2 + curl[(v m Vv hm ) - {v hm Vv m )\ . 
The term in the square brackets is equal to [{v m Vvh m ) — {vh m Vv m )] = [v, uj h ] and we obtain 

d t u h = [Vp, Vp]/p 2 + curl[v, u h ] . (71) 

For barotropic flows the first term in the r.h.s. becomes zero and we can see that o^/p is 
frozen-in field, 

/^\ = /^ \ (72) 



P J VP 

For H = u)h corresponds to the conventional hydrodynamic vorticity. 

In spite of the gauge dependence of the generalized vorticity, it frozenness gives us possibility 
to introduce the generalized helicity integral of motion. 

4.2 Generalized helicity. 



Now we can prove that generalized helicity, ha, defined in terms of the 'hydrodynamic' part 
of the velocity, 

h H = (u>h,Vh), (73) 

is integral of motion for barotropic flows. Namely, differentiating Eq. (|73|) and taking for account 
Eqs. (|70D , (fn]) we arrive for the barotropic flows to the local conservation law of the form (rather 
cumbersome calculations are presented in Appendix): 

d t h H + divq H = 0, = h H v + (x - v 2 /2)u h . (74) 

In analogy with the hydrodynamic case we arrive to the conclusion that the integral helicity 
Ih (defined by means of Eq. (|68|) ) is integral invariant, moving together with the fluid if the 
normal component of the vorticity tends zero, u>hn = 0, on the surface of the corresponding 
substantial volume V. Note that the condition Uhn = is invariant of the flow (due to the 
frozen-in character of i^h/p) and therefore it can be related to the initial surface only. 

4.3 Generalized Ertel invariant. 



Let us show here that there exists strict generalization of the Ertel invariant for the MHD 
case. For this purpose let us prove that without any restrictions related to the character of the 
flow the quantity 

h E = (w A , Vs) (75) 

satisfies conservation law of the form 

d t h E + divq E = 0, q E = h E v. (76) 
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Equivalently, the quantity = h E /p is transported by the fluid 

Da E = 0, a E = h E /p, (77) 

presenting a- type invariant. For the barotropic flows it immediately follows from the fact that 
u>h/p presents frozen- in field if one takes for account the composition rules given by Eqs. (|58|) and 
(|57D. In order to make the proof for the non barotropic flows more transparent let us consider 
something more general situation. Let J satisfy equation of motion of the form 

DJ = (J, V)v + Z, (78) 



differing from the frozen field equation flog ) by existence of the term Z that violates homogeneity. 
Then, if a represents any scalar Lagrange invariant, we have 

D(J, Va) = (DJ, Va) + (j, D(Va)) = (Z, Va) + (((J, V)v), Va) - (j, (Vv m ) ■ <9 m a)) . 
Here the two last terms cancel and we get 

D(J, Va) = (Z, Va) if DJ = (J, V)v + Z, and Da = 0. (79) 

For Z = this relations prove the generating rule of Eq. ([5*8]). But we can see that (J, Va) 
will present the local Lagrange invariant under more restrictive condition (Z, Va) = 0. That is 
the case for the Ertel invariant: Z = [Vp, Vp]/p 3 is orthogonal to Vs due to the fact that the 
scalar product of any three thermodynamic quantities is equal zero (because any thermodynamic 
variable in the equilibrium state can be presented as function of two basic variables). This ends 
the proof. 

The conserved integral quantity associated with a E is 

I E = [drh E , d t I E = 0. (80) 
Jv 

Note here that by the structure X E is not gauge invariant in contrast to the hydrodynamic case. 
Let us examine it change under gauge transformation changing =>- v^, vm =>- V M with 

Then 

l' E -l E = dr(Vs, J h - uj h ) = dr(Vs, u> M - oj' m ) . 
Jv Jv 

But (Vs, u>m — &m) — — div [Vs, {y' M — Vm)\ an d therefore we can proceed as follows 

X' E - X E = - f dE(n, [Vs, (v' M - v M )]) . 

JdV 

Taking into account that V M — Vm = — [h, M' — M] we obtain 

1 e -1e= [ dE(n,[Vs,[h,M'-M]]). (81) 

JdV 
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Inasmuch as both M' and M satisfy Eq. (11), their difference is governed by homogeneous 
equation 

d t M = cuW[v,M], 

i.e. m = M/p is frozen-in field. Then we arrive to the conclusion that the vector [Vs, [h, m]] 
entering the integrand presents frozen-in field, as it follows from recursive relations Eqs. fl5T|) - 
d59|). Therefore, if we adopt relation (n, [Vs, [h, m]])| a y. = as initial condition, then it holds 
true for all moments. But we cannot choose the (initial) substantial volume in such a way that 
relation 

(n, [Vs, [h,m']])|^ = (n, [Vs, [h,m]])|^ (82) 

holds true for any change of the gauge. Thus integral Ertel invariant is gauge dependent. 
Nevertheless, we can point out some subset of the gauge transformations under which is 
invariant. Namely, let M| t=to = /H for some initial moment, t = to, where / have to satisfy 
condition (H, V/) = 0, following from the divergence-free character of M. Then relation ( |52"D 
fulfills for the initial moment and therefore it holds true at all moments also. The specific choice 
/ = leads to additional restriction for the gauge transformations but it is convenient due to 
its simplicity. Summarizing, we can say that the Ertel invariant is partly gauge independent. 



14 



5 Conclusions. 



The results obtained can be summarized as follows. First, there is presented variant of 
introducing the canonical description of the MHD flows by means of the variational principle 
with constraints. It is shown that in order to describe general type MHD flows it is necessary 
to use in the generalized Clebsch type representation of the fluid velocity field vector Clebsch 
variables (the Lagrange markers and conjugate momenta) along with the entropy term (compare 
papers |l6j describing hydrodynamic case) and the conventional magnetic term introduced 
first in the paper | I7| . Such complete representation allows one to deal with general type MHD 



flows, including all type of breaks, see Ref. |14 |. Second, it is proved that introduced in the 



paper generalized Weber transformation leads to the velocity representation, which equivalent to 
the one introduced by means of the variational principle. Third, there is proved existence of the 
generalized Ertel invariant for MHD flows. Forth, there are generalized the vorticity and helicity 
invariants for the compressible barotropic MHD flows (first discussed for the incompressible case 
in cf. ||). Fifth, the relations between the local and integral invariants are discussed along with 
the gauge dependence of the latter. 

As a consequence of the completeness of the representation proposed we arrive to the correct 
limit transition from the MHD to conventional hydrodynamic flows. The results obtained allow 
one to deal with the complicated MHD problems by means of the Hamiltonian variables. The 
use of such approach was demonstrated for the specific case of incompressible flows in the series 
of papers |8], |J devoted to the nonlinear stability criteria. Emphasize, that existence of the 
additional invariants proved in our paper is of very importance for the stability problems. 

Note here that existing of the additional basic invariants of the motion makes it actual 
to examine the problem of the complete set of independent invariants, cf. Q. This problem 
needs special discussion together with related problem of their gauge invariance. One more 
open problem is connected with the great number of the generalized coordinates and momenta 
involved in the approach discussed. Here the question arises if it is possible to reduce this number 
without loosing the generality. 



Appendix A 



In order to prove Eq. (|32|) let us substitute J from Eq. (|31|) into expression [J, h] k dx k /dai. 
Then 

[J, h] fc ^ = [Dm, H] fc ^ - [(m, V)v, H] fe ^ = 
= ^D ([m, H] fc ) - ([m, D(ph)] h + [(m, V)v, H] fc ) ^* . 

Proceeding with the terms in the second brackets we obtain 

[m, D(ph)) k + [(m, V)v, H] fc = [m, h] fc ■ Dp + [pm, Dh] k + [(m, V)v, H] k = 
= -[M, h] k ■ divv + [M, (h, V)v] fe + [(M, V)v, h] k = -[M, h] s d k v s , 
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where M = pm and there are is taken for account dynamic equation Dh = (h, V)v and identity 

[M, (h, V)v] fe + [(M, V)v, h] fc = [M, h] k d s v s - [M, h} s d k v s . 
Introducing for brevity notation 

Y = [m,H] = [M,h], 
we can represent the r.h.s. of Eq. (jS3| ) as 

dx k v dx k^ dx k dv s dx k d ( dx k 

■ DY k + Y s —- d k v s = —— ■ DY k + Y s — = —— ■ DY k + Y s — (Dx s ) = D \ Y k 



oaj oa; ocl; Odj oa; oat. \ oa 



This proves Eq. (|32|). 

Let us check up now the integral relation (^). It is sufficient to prove the differential form, 
namely 

D(M,dE) = (J,dS), (85) 

where dH presents some infinitesimal oriented area moving with the fluid. It can be presented 
in the form 

d£= [dli,dl 2 ], (86) 

where dli, d\2 are frozen-in linear elements. Thus, dl a , a = 1, 2, are invariants of the J type and 
satisfy equations 

D(tfl a ) = (dl a ,V)v. 

Consequently, from the recursion relation Eq. ( pCf ) it follows that pdH is L- type invariant and 
therefore is governed by dynamic equation of the form: 

D(pdll) = -V(pvcffi) + [v,curl(pdE)], 

or in the coordinates, 

D(pdE i ) = -(pdE k )d i v k . (87) 
Now it is easy to prove relation ( |55D without any restrictions for the type of flow. Namely, 

D(M, dS) = £>(m, prfS) = (Dm, pdS) + miD(pdEi) = 
= (pdE, (m, V) v) + (J, tffl) - m iP dE k diV k = (J, dS) . ' ' 

In order to prove the helicity conservation, Eq. (0), let us consider some scalar quantity of 
the form 

Y = KJ), 



where J is frozen-in field. Then, taking for account that Eq. fl70D for the barotropic flows can 
be presented in the form 

Dv h = -V(x - v 2 /2) - v hm ■ Vv m , x = J dp/p, 

we obtain 

DY = (Dv h , J) + (v h , DJ) = - (V(x - v 2 /2), J) . 
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For J = uj h / p we proceed 

D(v h , w fc /p) = -p- 1 (V(x - ^ 2 /2), w h ) = -p- l div ({x - ^ 2 /2K) • 

Then 

D(v h , oj h ) = pD(v h , w h /p) + (v h , uj h /p)Dp = -div ((% - v 2 / 2 )^) - (v h , u h ) divv, 

or 

d t (v h ,u) h ) = -divqh, q h = (x - v 2 /2)u: h + v (v h , u h ) (89) 

that evidently coincides with Eq. (]7^) . 

It is noteworthy that the proof is valid for any J- type invariant if the field pJ is divergence- 
free: 

d t {pJ, v h ) = -dzuq , q = (x - w 2 / 2 )P J + v (P J > v /i) for div(pJ) = 0. (90) 

For instance, choosing J = h immediately leads to cross-helicity invariant if one takes for account 
that (H,v fc ) = (H,v). 
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